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<D ■ Abstract 

m 

A field-theoretical model for non-singular global cosmic strings is presented. The model 
is a non-linear sigma model with a potential term for a self-gravitating complex scalar field. 
Non-singular stationary solutions with angular momentum and possibly linear momentum are 
obtained by assuming an oscillatory dependence of the scalar field on t, ip and z. This dependence 
has an effect similar to gauging the global (7(1) symmetry of the model, which is actually a 
Kaluza-Klein reduction from four to three spacetime dimensions. The method of analysis can be 
^/-^ ' regarded as an extension of the gravito-electromagnetism formalism beyond the weak field limit. 

f"^ , Some D — 3 self-dual solutions are also discussed. 

^ ' PACS: 11.27.+d, U.lOLm, 04-20. Jb, U.lOKk, 98.80. Cq 

O ' 

1 Introduction 

Ph. 

Cosmic strings are a special kind of topological defects [I] , which could have been formed in a symmetry 
breaking phase transition in the early universe. Often, the broken symmetry is assumed to be local, 
but the possibility of global strings has also been extensively discussed pQ. The prototype is Cohen and 
Kaplan's solution E3 E] of a complex scalar with a U(l) symmetry breaking potential. The string 
solution turns out to have a repulsive gravitational field with a curvature singularity at a finite distance 
from the axis. The singularity is unavoidable as was shown by Gibbons et al. 0] for the standard 
Lagrangian of a complex field with a non-negative potential term, under the assumptions that the 
solutions are static and have symmetry under boosts along the string axis. This result may be regarded 
as the analog in the gravitational case of the D = 2 + 1 dimensional version of Derrick's theorem |HJ, 
excluding the existence of static global strings with finite energy per unit length in flat space. 

It was shown explicitly that the singularity of the Cohen-Kaplan solution can be avoided if time 
dependence is allowed, which results a string solution with de-Sitter like expansion along the string axis 
[Hj. An attempt to get a static non-singular global string was tried [J] by relaxing the symmetry under 
boosts parallel to the string axis (z), but without relaxing the corresponding condition (T t * = T*) in the 
string core, so the solution found there has lost its interpretation as a global string gravitational field. 

There are additional ways to overcome the "no-go" theorems and find string solutions. All of them 
involve modifications of the Cohen-Kaplan model. Probably the most well-known is gauging the U(l) 
symmetry group to become local, thus obtaining the Abelian Higgs model. The main consequence 
of this modification is that the long range massless Goldstone field transmutes into the longitudinal 
component of a short range massive vector field. This model admits string solutions (flux tubes) in flat 
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space and in the self-gravitating case. However, we will not discuss this possibility here since almost all 
that can be said about it, has been already said Other ways are either restricted to the flat space 
case, or to the gravitating case or are valid in both. One is the Q-ball solution [S], which overcomes 
Derrick's theorem by allowing oscillatory time dependence of the complex scalar, so as to avoid the 
staticity assumption without introducing time dependence into the energy-momentum tensor and other 
observables. In addition, the potential term is modified in such a way that it has another symmetric local 
minimum. From a more physical point of view, these modifications result a non-vanishing value of the 
conserved global U(l) charge, which stabilizes the spherical solutions against collapse. The string-like 
analogue ("Q-string") is similarly stabilized by a finite global U(l) charge per unit length (see however 
[Hj for a different situation). Spinning versions have been also recently considered ^Oj- Adding gravity 
into the Q-ball system does not change the overall picture of the solutions, but rather yields the so-called 
Q-stars [TT] . 

A second way which is also based on the existence of a global U(l) charge but in which gravity is 
indispensable, gives the boson stars [T^llT^lHl) . 

It is also possible to get similar "non-topological solitons" in systems with more than one complex 
scalar field |14j , but since we want to stick to a minimal field content in our model, we will not discuss 
those as well. 

The simplest way to add a topological charge is to assume for the scalar field a "hedgehog config- 
uration" , where it tends asymptotically to symmetry breaking minima of the potential, mapping this 
way spatial infinity into the vacuum manifold. In the original U(l) system, static string-like solutions 
fall indeed into homotopy classes characterized by an integer winding number, but such solutions are 
excluded by the "no go theorems" mentioned above. However, we may think of a second version of a 
"hedgehog configuration" , where we turn the target space from a trivial vector space into a compact 
manifold of a non-linear sigma model. Since we are after string-like solutions, the scalar field defines a 
mapping to target space from the compactified R 2 , which is transverse with respect to the string axis. If 
we take the simplest target space S 2 , the kinetic part of the Lagrangian will be the usual one of the 0(3) 
non-linear sigma model. The potential term will necessarily break explicitly the 0(3) symmetry, but we 
will assume it is still {7(l)-symmetric. It will be also taken until further notice to be non-negative with 
minima at U = 0. In this form of hedgehog configuration, the scalar field should tend asymptotically 
to a J7(l)-symmetric minimum in order to correspond to a linear defect, which will be both localized 
enough and have a topological charge. 

To make the discussion more concrete, we write down the action in the general form 

S = J d D x\/\g\ (|*|)V M **V<I> - U(\<P\) + i-fl) (1.1) 

where n is the D dimensional gravitational constant (we will discuss only D = 3,4) and £(|$|) is a 
non-negative dimensionless function, which serves as a Weyl factor of a conformally-flat target space 
metric of the non-linear sigma model. The S 2 target space is obtained by 

g(l$l)= (i + IW) 2 (L2) 

where (i is the diameter of the 2-sphere and sets the energy scale. In D dimensions it has dimensions of 
energy to the power (D — 2)/2. More conventions: V M is a covariant derivative, signature (+, —,...,—) 

and the Riemann tensor is J?xuw = ®v^\u ~ ^^Xv + 

Note that the "linear" Higgs (or rather Goldstone — but we will use the former term) system is 
recovered either by taking £ (|$|) = 1, or using the sigma model Lagrangian in the limit of small fields 
or formally for fi — > oo. 

In flat spacetime this model has well-known static string-like solutions, but only with a vanishing 
potential - again a simple consequence of Derrick's theorem. This is the original 0(3) sigma model 
We will still use the term "sigma model" for the action (| 1 . 1 p . but refrain from characterizing it by the 
symmetry group, since it is (explicitly) broken by the potential term to U(l). 

The flat space string-like solutions were found to be unstable and may either shrink or spread out 
|151 116). They may be however stabilized by allowing (as above) an oscillatory time dependence and 
adding a potential term. This produces the so-called "Q-lumps" , which are stable rotating topological 
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string- like solutions ^| E|. As described, these solutions have oscillatory time as well as angular 
dependence (around the string axis): $ = ^,f{r)e lljnip ~ ult \ where m is an integer and uj real. The 
non-vanishing global U(X) charge results from the oscillatory time dependence, while it is the mixed 
angular and time dependence which produces the non- vanishing angular momentum density, T* The 
topological winding number is m. 

The 0(3) sigma model may be also stabilized by a gravitational field in the form of a cosmological 
expansion J^j. Closer to our present purpose is the Comtet-Gibbons solution [201, which is a static 
cylindrically symmetric (self-dual) solution of a self-gravitating 0(3) sigma model (with a vanishing 
potential). 

In this paper, we fill some of the gaps which are evident from the above discussion, and find new 
non-singular global string solutions, which are topologically non trivial: self-gravitating Q-lumps and 
generalizations which introduce also a momentum current along the string axis - a possibility which 
has not been considered before. This is reflected by T' ^ in case of rotation or T' ^ in case 
of translational motion along the string axis. This will be done by taking the more general oscillatory 
behavior $ = nf{r)e l ^ qz ^ wtJrm ' p \ Some non-topological solutions like self-gravitating spinning Q-strings 
will be discussed in a future publication. 

The expectation that one may get non-singular stationary gravitational fields in this way, is based on 
the observation that introducing an oscillatory z and/or i-dependence is consistent with the existence 
of two commuting Killing vectors, which generate translations in time and along the string axis, and 
has the effect of gauging the global symmetry and turning the global sigma model into a gauged one. 
This way, the original D — 4 gravitating global sigma model is reduced in a Kaluza-Klein way into a 
D = 3 gauged sigma model coupled with D = 3 gravity and electromagnetism (and a dilaton) . Since it 
is known that the gravitating Nielsen-Olesen flux tube produces a non-singular gravitational field, even 
when a dilaton coupling is allowed 1221 |2~3"| , there is a good chance it will share this characteristic 
with the type of solutions described above. As the dilaton introduced by the Kaluza-Klein reduction is 
also a long range field, the gravitational field of the global string solutions obtained in our model will 
be finite but not asymptotically flat. Only the pure D = 2 + 1 system discussed in sec. Elwill admit 
asymptotically flat global vortices. 

We proceed and end this section by writing the field equations of our model. They are the following 
generalization of Klein- Gordon equation 

emw* + + ^ = o (i.3) 

together with the Einstein equations (in term of the Einstein tensor G^) 

- —Gfj, u = v = [d^*d^ + d^*d^ - (V A $*)(V A $).g^] + 17(1*1)^ (1.4) 

which are equivalent to 

\r*» = -^um)^ - id^*d„<P + wa^] (i.5) 

2 Field Equations for Sigma Model Global Strings 

We start with the Z?-dimensional action (|l.l|l and look for string-like solutions, where the metric tensor 
is written in the form 

ds 2 = g^{x x )dx»dx v = e 2a{xl) h a p{x' 1 )dx a dx p - e^-^^^dz + A a {x' 1 )dx a ] 2 . (2.1) 

The indices a, j3... go from to D — 2, and we denote x D ~ l = z. This is just the Kaluza-Klein reduction 
from D to D — 1 dimensions. If we also take $ = ip(x a )e lqz , the Z)-dimcnsional gravitating sigma model 
with a global U(l) symmetry turns into a gravitating gauged sigma-Maxwell-dilaton system in D — 1 
dimensions, as it is evident from the Lagrangian which becomes: 
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h\ (~£ m)[D a rD a ^ - q 2 e 2( * D ~ 2 ^\ 2 ] - e 2 °U{\^ + 
(R{h) + {D- 2){D - 3)V Q aV Q (7 - i e ' 2{D ~ 2)o ' F aP F aP + 2(D - 2)V Q V Q ^ (2.2) 

after we define a usual covariant derivative D a = d a — iqA a and field strength F a p = d a Aj3 — dpA a . 
The last term can obviously be omitted. All geometric quantities are related to the (D — 1) metric h a p. 
In the next step we assume stationarity, i.e. parametrize the scalar field and the D — 1 metric as: 

ip = <t>{x k )e~ lu)t , h a p{x 1 )dx a dx f} = N 2 (x k )[dt + L^x^dx 1 ] 2 - lio {x k )dx l dx j (2.3) 

where the indices go from 1 to D — 2. Now we do not write the Lagrangian or the action in the 
D — 2 dimensional notation, but assume D = 4 and write directly the field equations in terms of the 
following: 

• Vi the covariant derivative with respect to the 2-dimensional metric 7ij , 

• Rijij) the corresponding Ricci tensor, 

• Lij = diLj - djLi = ly/\i\eij. 

• F oi =NE i , F 0i = -D i /N 

. Fa = -By/Whj, f* = -zretf/Vffl 

• Dj = dj — i[qAj — (ui + qA n )Lj] = dj — iqAj 

. J fc = -^(|0|)(0*^-05^*) 

We use the usual convention of raising and lowering indices by jij and its inverse. Note the distinction 
between D l which is an "electric" (displacement) field (obeying Di = E^ + N H \/\j\eijlJ), and D l 
which is a covariant derivative. We will also use E 2 to denote EiE 1 . This way of dealing with the field 
equations is essentially an extension of the gravito-electromagnetism formalism beyond the weak field 
limit, where it is usually used |241 125j . The components of F a p will be referred to as gravito-electric 
and gravito-magnetic. 

Einstein equations split therefore according to the components (00), (zz), (iO), (iz), (Oz) and (ij). 
The first 5 are: 



1 tt2 "p 

_ Vi (JVVV) + — ' 



V i (iVVV) + e~ 4CT = K{e 2 °Um+q 2 e^\<t>\ 2 £{\4>\)) (2.5) 



c 



dj(N 3 £)-- ^=—E,j = 2n{u + qA a )J l (2.6) 



djiNHe'^) = -2nqJ l (2.7) 



1 T7 f„-4vr>i\ _ n.,„(r fi " + qA Q ^ 



-V,(e- 4 ^) = 2 Kq lL i r-— ] ^\cl>\'£(\ct>\)) (2.8) 



and the (ij) components are: 



1 N 2 £ 2 1 

%(7) + ^ViV,-iV - —7« + MiadjO + -(E 2 7y - - E^e 



-k ((2e 2 °U(\4>\) + q 2 e 4 °\cj)\ 2 £(\cf>\)) 1%3 + ^SQ^D^D^ + D^D^ (2.9) 
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The equation for the scalar field is 

Dm + ' * A#OV + (^g^ - ^) ( i + ♦ - £ « , . o (2.10) 



Another equation which is not independent, but will be useful is the (ij) component of (|1.4|) 
1 N 2 £ 2 

jrCViVjN - 7 tf V*V*JV) — 7ij + 2<V<V - ^ M d k adi a + 

4,7 / H 2 + F, 2 



2 V 2 

2 1 



| (f^DS^V + 2e 2 °UM) + (q 2 e*° - { ^pL \ \^S(\cf>\^ 7ij (2.11) 

3 Three-Dimensional System 

In order to get a clearer picture of the system, it is useful to deal with a lower dimensional case if 
possible. Although this set of field equations is genuinely four-dimensional, it contains also in a special 
limit the field equations for the D = 3 system. These are equations 1)2. 4p . (|2.6|) . (12.91) . (|2.1U|) with the 
substitutions: q — 0, a — and ^4 a = 0, which imply J3j = Di = and B = H = 0. We get therefore: 

1 /v 2 / 2 / ^ 2 \ 

±ViVW+— = -^2[/(|0|)-— |0| 2 £(|0|)j (3.1) 



&(A^) = 2^/' (3.2) 

1 N 2 ^ 2 / 1 \ 

%(7) + ^V i V i iV-^-7y = - K ^2[/(|0|) 7y - + -SMXDiPDjt + DsPDnt,)) (3.3) 

and for the scalar field: 

Note that a solution of the D = 3 system is not automatically a solution to the D = 4 system with the 
above-mentioned substitutions. The main reason is eq l|2. 5J1 which gives an additional condition, which 
is consistent with the D = 3 theory only in the absence of self-interaction (i.e. U — 0). On the other 
hand, the fact that there is an explicit dependence in the Einstein equations upon the dimensionality 
D, does not hinder the possibility to use them for the D = 3 field equations. Accordingly, the D = 3 
analogue of eq (|2.11() is: 

1 N 2 f 2 K 

-{ViVjN - 7y -V fe V fe A0 - — 7y = --SQtMDiFDrf + D^D^) + 

| ^(101)^0*^0 + 2U{\4\) - ^rl0| 2 £ (101)) 7tf (3-5) 

4 Self-Dual Solutions 

It is known that the field equations simplify considerably if self-duality is possible (see e.g. ref. I2H1 and 
references therein) . So we will start in the simplest case of self-dual solutions in the D = 3 system. The 
self-dual solutions obey 

N = l, D i <t> = ir)-^=D j <f> (4.1) 

vm 
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where r] = ±1 corresponds to self-dual or anti self-dual solutions. In the following we will refer to both 
as "self-dual". This is consistent with eq. (|3.1|) or (|3.3|) . provided the following relation holds: 



— + 2UM)-LO 2 \<j)\ 2 £(\cj ) \)=0. 
Zk 



(4.2) 



This relation is also consistent with eq (|3.4|l . Since the current also simplifies considerably for self-dual 
solutions, we get from eq 13.211 

d i £ = 2r l KUj\<t>\£(\<j>\)d i \<p\ (4.3) 

This equation can be integrated to give £(\4>\) for any given £ (\4>\). In the most familiar case of the S 2 
sigma model defined by eq. (|1.2I) , the integration is trivial and we find the following form of potential 
term, which will give rise to localized self-dual solutions: 



UsdM) 



^V(|$| 2 //i 2 ~^ 2 /2) 
2(1 + |$| 2 //^ 2 ) 2 



(4.4) 



We chose the integration constant such that £ will vanish asymptotically i.e. for fields that satisfy 
$| — > oo as r — > oo. It is interesting to note that the linear sigma model (£ (\4>\) — 1) needs for 
self-duality a potential which is unbounded from below, which is a property shared by other systems 

129 - 

Since target space is a sphere, it is natural and useful to use an "angular" field defined by |$| = 
/itan(0/2). In terms of this field, we get for example £ — cos 4 (9/2) and see that |$| — > oo is actually 
just the "south pole" of target space, = tt. Thus, the self-duality potential can be rewritten as: 



c/ SD (e) = ^(sin 2 e 



2kh 2 cos 4 (0/2)) 



(4.5) 



This generalizes the flat space self-dual Q-lumps to be self- gravitating. The flat space limit is 

reproduced by taking k — 0. From a different direction, this can be regarded as a spinning generalization 
of the self-dual Comtet-Gibbons solutions |2()j . Those are reproduced for oj = 0. Explicit solutions are 
presented in sectional 



5 Rotational Symmetry 

We now return to D — 4 and specialize to rotational symmetry. We are left with the 2-metric 

<l.r- = a 2 {r)dr 2 + (3 2 {r)dip 2 (5.1) 

where we postpone the gauge fixing to a later stage. The vector A a will have only two non-vanishing 
components Ao(r) and A v (r), from which we derive the gravito-electric and gravito-magnetic fields 
E r (r) and B{r). Similarly, Li will have only one non- vanishing component namely L v (r); thus £ will 
also depend on r only. We will also use A v defined by qA v = qA v — [w + qA )L v . 

We define a dimensionless potential function u by U = A/Lt 4 u(0) where A is a dimensionless coupling 
constant (in D = 4), and a dimensionless Newton's constant 7 = Kfi 2 and get therefore the following 
set of field equations for the S 2 model. From equations H2.4|) - (|2.8|) we get 



1 



//37V'Y iV 2 ^ 2 



H 2 



Na(3 V a 
1 fPNa'V H 2 -{E r /af 



Na/3 V a 



-4a 



-4a 



Na 



{N 3 £)' 



2 

0N 



(m + qA Q f 
N 2 
,sin 2 9 



g 2 e 4CT 



sin 2 6 



2Xn 2 e 2a u{Q) 



2>^ + A ^\ ( e) 



q e 



HE r e 



1 



Na/3 \a 



Na 



(NHe 



Hte 



--(ui + qA )(m - qAtp) sin 
-^q(m — qA v ) sin 2 



7 
'27V 2 



q{uj + qA ) sin 2 



(5.2) 
(5.3) 

(5.4) 
(5.5) 
(5.6) 



G 



where the last one is a combination of (|2.7(l and 1)2.8(1 . From eq 1)2. 9JI we take the equation for i?(7) 



2 //3'Y 3A^ 2 H 2 + E 2 /a 2 _ ia a' 2 



a(3 \ a J 2 2 a 2 

/ 6^ / (m-gi y ) 2 (c + gAo) 2 2 4 \ sin^e 



-7 ^ + ^^F^ + + S^j — + j (5.7) 

and the (tpip) equation: 

1 W NH2 , V^ff (m -]^ + ^^ + 2V^(e)) (5.8) 



a/3iV V « / 2 2a 2 ' VV /3 2 

For the scalar field we get: 

1 fN(30'\' f{uj + qA a ) 2 (m~qA v ) 2 2 4 \ sin(26) 2 2(T du 



and we have also to add the following relations: 

* = -#. ^=^(^o-A;), £=|. (5.10) 

For completeness we also note that B = —A' v /af3. 

Note that once we make a gauge choice, one of these equations becomes dependent on the others. 
When we solve them numerically we use eq. (|5.8|) as a check. 

Few soliton-like solutions to this system are already known in the gravitating case - all of them 
static (u> = q = 0, L v = A v = Aq = 0). For the "linear" case (£(|$|) = 1) with a potential with a 
U(l) symmetry breaking minimum, there exists Cohen and Kaplan's global string 2 which is singular as 
already mentioned in the introduction. In the non-linear sector, the Comtet-Gibbons string-like solution 
PU] is a self-dual solution with vanishing potential. 

Remotely related solutions are the dilatonic-Melvin solution [271 I22| . and the three-dimensional 
dilaton-Maxwell-Einstein system that were discussed by some authors [281 1291 13UI 1ST] . The main dif- 
ference between our model and the above mentioned dilaton-Maxwell-Einstein papers, is that none of 
them has suggested a possible source for the angular momentum or electromagnetic fields of those solu- 
tions. It may be thought that those dilaton-Maxwell-Einstein papers may be considered relevant to our 
model in the asymptotic regions, where the scalar field settles to its vacuum value, but even from this 
point of view there are significant differences. The dilaton in the present system is a metric component 
and not an additional field, so its coupling to other fields is not arbitrary. Similarly, the "electric" or 
"magnetic" fields are also of gravitational origin, and since they have specific source, the field equations 
are inconsistent with the separate conditions E r — or H — 0. Consequently, the solutions presented 
here cannot be classified as either electric or magnetic but always (for non-vanishing self-interaction) 
have both kinds of fields or none of them. A third difference is that most of the above-mentioned papers 
introduce a non-vanishing cosmological constant, which we avoid here. 

We solve the equations numerically taking the "Kasner gauge" a(r)e cr<r -' = 1 so that g rr = — 1. 
In order to get localized solutions, we solve the equations with the boundary conditions which ensure 
regularity at the origin and will be as close as possible to asymptotically vacuum. That is, we impose 
on the 3-metric 

N(0) = 1, N'(0) = 0, 0(0) = 0, P'(0) = 1, M°) = °> ^(°) = - ( 5 - U ) 
As for the gravito-electromagnetic components, ^4o(0) can be any finite constant while 

4,(0) = 0, lim A (r) = 0, ^(0) - (5.12) 

r — >oo 

(for ^(oo) see below) and for the "dilaton" and the scalar field 

cr(0) = 0, cr'(0) = 0, 6(0) = 0, lim 6(r) = 9 (5.13) 
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where 9o is the location of the minimum of the potential. The value of 9o distinguishes between two 
types of hedgehog configurations: a "sigma model hedgehog" and a "Higgs hedgehog" . The first can 
be realized if the scalar tends asymptotically to the south pole, i.e. Qq = ir. In this case ^(oo) can 
be any finite constant and the gravito-magnetic flux is not quantized. The second corresponds to any 
other value < 9o < tt which is associated with a quantized flux: A v (oo) — m/q. In order to discuss 
both kinds of solutions, we may use the following simple potential function: 

/ sin 2 (9/2)-sin 2 (9 /2) 
Me)= { SnW^ 

where 9o is, as mentioned above, the location of the potential minimum. It is actually the standard 
Higgs potential multiplied by the conformal factor ^d^Q and written in terms of the angular field. The 
case Qq = tt corresponds to a "single well potential" . We will also use the "double well potential" , which 
is a special case of 1)4.51) with two minima at the two poles 9 = 0, tt: 

M 2 (e) = sin 2 e (5.15) 

These minima do not break the symmetry so they allow topological solutions of the "sigma model 
hedgehog" type. This is the potential used by Leese in D = 3 flat spacetime for obtaining the 
self-dual Q-lumps. Ward ^S] used another form given by u(9) = sin 2 0(1 + cos 2 9), which we will not 
use here. All the other potentials are shown in figure 



(5.14) 



6 Special Solutions: Strings with Linear and/or Angular Mo- 
mentum 

We solve numerically the field equations, using the MATHEMATICA package, after rewriting them as a 
system of 14 coupled first order ordinary differential equations. We use the dimensionless radial variable 
x = //r, and rescale accordingly the other dimensional quantities q//i — > q, to/ \x — » oj, nf3(r) — > (3{r) and 

The following sections describe the main characteristics of the solutions. In all cases we take m = 1. 



6.1 Higgs-like Potentials 

We start by discussing a family of solutions for potentials of a Higgs type, whose minima break spon- 
taneously the 17(1) symmetry - eq 1)5.14(1 . The main disadvantage of the Higgs type models is that the 
solutions are not localized "enough", so the gravitational field is still singular, and we present them here 
only as a basis for improvement. Figure [21 depicts the fields of a static solution for the case 0o = ir/2. 
It is very similar to the Cohen-Kaplan solution. For small 7 the singularity appears for very large r (or 
x - for 7 = 0.2 and u> = the singularity is already at x ~ 10 9 ), so for demonstrative reasons we took 
quite a large 7. There are two kinds of singularities. The first is the one appearing in figure [21 where 
N, e~ a and j3 vanish while fie" diverges for some finite r. For u> — Ne a — e~ a so N and e~ ff vanish 
at the same point. If \ui\ > but smaller than a certain bound (which depends on 7), Ne a decreases 
faster than e~ CT , so it is the zero of Ne a or rather N which determines the singularity. If however \lo\ is 
large enough, (5e a changes its trend and vanishes while N and Ne a diverge. This last behavior is seen 
in figure [21 A similar behavior is found for q ^ as well. 

This behavior is a direct result from the existence of the long range Goldstone fields that exist in 
this model. This is the origin for the singularity of the static metric of the global string, and the other 
long range effects in other global defects. 

As in the Cohen-Kaplan solution, it is the contribution to the energy density from the angular 
dependence which does not vanish fast enough to produce a localized source. In addition to the angular 
dependence, we allow here also time dependence and z-dependence which yield similar contributions all 
present in the right-hand-side of H5.2Jl - l|5.8|l . Since all these terms are proportional to sin 2 0, it is quite 
easy to identify the way to improve the asymptotic behavior namely, taking 0o = tt or |$| — ► oo as 
r — > 00. Note that the opposite direction, O << 1 gives the Higgs limit, as already mentioned. 
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6.2 Single Well Potential 

In order to have solutions with Oo = ir (or |$| — > oo), we take the following simple monotonically 
decreasing potential 1/(1$!) = A/i 4 /(l + $| 2 //z 2 ), which is one of the "vacuumless" potentials of Cho 
and Vilenkin (22123121) who also found only singular string-like solutions. However, since we use it 
within the sigma model framework, this potential, unlike those of Cho and Vilenkin, is not vacuumless 
but a single-well potential with a single minimum at the south pole (O = n) as is evident from its 
behavior in terms of the angular field: 

u(Q) = cos 2 (0/2) (6.1) 

Thus, we may get non-singular global strings and indeed we find non-singular stationary solutions, 
though not static ones. There are two kinds of string-like solutions: purely rotating (PR in what 
follows) where no gravito-electromagnetic fields exist, and rotating solutions with additional gravito- 
electromagnetic fields (RGEM). On a first sight, one may think that the fields of the latter may be either 
gravito-magnetic or gravito-electric, according to the field components produced by the momentum 
along the string axis. However, this is not the case (if self-interaction exists) since the field equations 
are inconsistent with either E r = or H = 0. Therefore, only solutions with gravito-electromagnetic 
fields exist. 

Both kinds of solutions (PR and RGEM) have w / and the PR solutions have q = 0. On the 
other hand, RGEM solutions exist for either vanishing or non-vanishing values of q. The fact that q — 
allows non-vanishing gravito-electromagnetic fields is a direct result of the "Maxwell" equation for the 
"displacement" field D l eq (|2.8(l . which shows that it is D l (rather than E l ) that vanishes if q = 0. 
Since there is no proportionality between these two gravito-electric fields, = does not imply E % = 
nor B = (or H = 0) but only a specific relation between them, somewhat similarly to a medium 
with electric and magnetic permeabilities. Therefore, these fields can be non-zero even when q = 0. 
Another way to see that there should exist RGEM solutions when q = is from the fact that it is 
possible to "gauge away" q by a coordinate transformation which however will not transform to zero 
the gravito-electromagnetic fields 1 . The case q = splits therefore into two branches: a PR branch and 
a RGEM branch. 

The fields of a purely rotating solution are depicted in figure 0J There is only a narrow range of 
u values, outside of which only singular solutions exist. Those of a rotating gravito-electromagnetic 
solution are presented in the two parts of figure RGEM solutions with q = are similar and are not 
shown. The fact that RGEM non-singular solutions exist with smaller uj is a general one, and may be 
understood by the fact that less angular momentum is required to balance the scalar and gravitational 
self-attraction, in the presence of gravito-electromagnetic fields. Another significant difference between 
the two kinds of solutions, is the asymptotic behavior of e~ ff . While it is monotonically decreasing 
function in the PR case, it is monotonically increasing for RGEM solutions. This difference is a direct 
consequence from the field equations: it is obvious from equation l|5.3(l that for vanishing gravito- 
electromagnetic fields, the function <j(r) is monotonically increasing so e~ a ^ is decreasing. On the 
other hand, for RGEM solutions, the magnetic contribution is the dominant term for large distances 
and it comes with an opposite sign, so e~ CT ( r ) has the opposite large distance behavior - although it may 
start decreasing near the axis. 

Another kind of potential which we can consider is 

u(0) = cos 4 (6/2) (6.2) 

which is just the special case ©o = tt of l)5.14|l . We find non-singular PR solutions, which are quite 
similar to those of the previous single well potential, thus we do not present them. However, we were 
able to obtain only singular RGEM solutions for this potential. 

6.3 Double Well Potential 

We take the double well potential (|5.15|l . The two minima do not break the symmetry, so they allow 
topological solutions of the "sigma model hedgehog" type. The fields of a purely rotating solution appear 

1 we thank the referee for pointing out towards this direction. 
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in figure El As for the previous potential, there is a range of uj values outside of which only singular 
solutions exist, while 7 can be quite large. The fields of a rotating gravito-electromagnetic global string 
are presented in figured As before, the q = RGEM solutions are similar as they are related to the 
q ones by a coordinate transformation. The gravito-electromagnetic solutions become closed above 
a certain critical value of 7. The behavior of e~ cr is also analogous to that of the single well potential: 
decreasing for the PR case and increasing for the RGEM case. 

6.4 The Self-Dual Potential in D = 3, 4 

Now we take the potential (|4.5|) which, as we saw already, allows D — 3 self-dual rotating strings or 
rather vortices. We write it as in all other cases in terms of the dimensionless potential function it(0), 
but with the D = 3 coefficient 2 as U = A/i 6 u(9). The necessary condition for self-duality is therefore 
oj 2 = 8Xfi 4 . In order to obtain vortex solutions, we may directly solve the second order system of 
equations l|5.2jl . I|5.4jl . (j5.7jl and (j5.9jl with the substitutions q — 0, N = 1, a = 0, A = and A v = 
which imply E r = and H — 0. Alternatively, we may solve a first order system which include (among 
other equations which we do not present), the equations l|4.1[) and l|4.3|l used for rotational symmetry. 
A typical solution is shown in figure |SJ 

Now we can proceed and search for solutions with the same potential (up to powers of fi) in D = 4. 
It is not surprising that we don't find self-dual solutions. Yet there exist non-singular solutions, which 
we do not show since they are quite similar to those of the double well potential. 

6.5 No Self-Interaction 

For 17 = we reproduce the well-known Comtet-Gibbons string- like solutions [201) which have asymp- 
totically a conical geometry. Figure [5] shows the fields of a solution with a very large value of 7 that was 
taken in order to stress the angular deficit, which is directly seen from the asymptotic slope of j3(r). 

A new type of (purely) rotating solutions is obtained for non- vanishing u>. The fields of one of these 
solutions are depicted in figure ITU1 All these solutions have asymptotic cylindrical two-geometry, while 
the metric component N(r) increases linearly as r — > 00. Solutions exist for any \uj\ > 0. The limit 
lj — > is singular, and does not reproduce the conical uj = solutions. 

The fields of a rotating gravito-electromagnetic global string are presented in the two parts of figure 
1111 A non-rotating solution with gravito-magnetic field is also found and presented in figure IT^l Because 
uj = we have here, as in the Comtet-Gibbons solution, N — e~ a but now different from 1. This kind of 
purely gravito-magnetic solution is special to the case of vanishing potential and cannot exist otherwise. 

7 Conclusion and Outlook 

We have shown that non-singular global cosmic string solutions exist in a simple sigma model, if we 
also allow for rotation or momentum current. Another result is the self- gravitating generalization of the 
D = 3 self-dual Q-lumps. Further investigation is required in order to fully understand these solutions 
and their consequences. The geometry of the solutions should be clarified by geodesic analysis, and 
their possible role in a cosmological context also deserves further study. Another direction is a higher 
dimensional generalization of these solutions to be codimension 2 branes. These may be relevant to the 
extra dimensional resolution of the hierarchy problem, which was advanced recently by several authors 
OH! El OH] using similar kind of solutions. 

2 Note that ny? = 7 is still dimensionless as in D = 4 (or higher). 
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Figure 1: Dimensionless potential functions: a "Higgs-like" Ui(6) with O = t/2 - dashed, the "single- 
well" ui (6) with 6q = it, the double well «2(0) and the self-dual usr>(6). 
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Figure 2: The fields of a static global string solution for a Higgs-like potential. 
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Figure 3: The fields of a purely rotating global string solution for a Higgs-likc potential. 
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Figure 4: The fields of a purely rotating global string solution for a single well potential. 
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Figure 5: The fields of a rotating gravito-electromagnetic global string solution for a single well potential. 
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Figure 7: The fields of a rotating gravito-electromagnetic global string solution for a double well poten- 
tial. 
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Figure 8: The fields of a rotating self-dual global string solution in D = 3. Note that N(x) = 1. 
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Figure 9: The fields of the Comtet-Gibbons (static) global string solution. 
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Figure 11: The fields of a rotating gravito-electromagnetic global string solution without self- interaction. 
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